Directional emission from weakly eccentric resonators 
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It is shown that when a circular resonator is deformed in a nonintegrable way, a symmetry 
breaking of escaping rays occurs which can dramatically modulate the outgoing wave even for 
small perturbations. The underlying mechanism does not occur in integrable models for which 
the ray families can be computed exactly and is described in this Letter on the basis of canonical 
perturbation theory. Emission from deformed resonators is currently of immense practical interest in 
the context of whispering-gallery optical resonances of dielectric cavities and the approach outlined 
here promises simple analytical characterisations in the important case of small deformations. 
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Short-wavelength approximations have proved invalu- 
able in understanding and predicting the properties of op- 
tical resonators [1, 2]. They have in particular provided 
essential insight into the directional emission patterns 
that are observed when resonators are sufficiently asym- 
metric that ray dynamics display chaotic as well as reg- 
ular behaviour [3] , where features of wave chaos such as 
chaos-assisted tunnelling [4-6], scarring [7, 8] and dynam- 
ical localisation [6, 9] have been shown to play a role. Sur- 
prisingly, emission characteristics in the apparently be- 
nign limit of very weakly asymmetric resonators cannot 
readily be described using this approach [10]. The reason 
is that emission is essentially a tunnelling phenomenon 
whose short-wavelength approximation demands that we 
extend the underlying ray families into complex space. It 
has been established, however, that natural boundaries 
generically intervene [11], even in very slightly perturbed 
problems, which prevent analytic continuation as far as 
needed in the complex domain [12, 13]. In other words, 
the ray-dynamical data demanded by short-wavelength 
approximations of emission simply does not exist. 

In view of the current technological interest in opti- 
cal resonators [2], simple analytical approaches to this 
problem are clearly of great value. We outline one such 
solution in this Letter which works by substituting for 
the nonexistent classical data approximate solutions that 
can be extended as far as required in the complex do- 
main while remaining sufficiently accurate for use in 
WKB approximations. A related approach has recently 
proved successful in predicting tunnel splittings in near- 
integrable potentials [14]. Although restricted to rela- 
tively small deformations, the solution is capable of pre- 
dicting strongly directional emission patterns. 

A feature of the solution is that there is a dramatic 
difference between the behaviours of integrable and non- 
integrable deformations. This is true even of problems 
whose Poincare plots look quite similar. To illustrate 
this, emission patterns from two perturbations of a circu- 
lar resonator, one integrable and the other nonintegrable, 
are compared in Fig. 1. The details of these systems 
are described later but here we simply point out that 
even though the Poincare plot of the nonintegrable prob- 
lem seems to deviate much less from the perfectly cir- 



cular limit, it has a dramatically more directional emis- 
sion pattern. The enormous qualitative difference be- 
tween integrable and nonintegrable systems, irrespective 
of any gross similarity of the real ray dynamics, indicates 
that simple geometrical characteristics of the deforma- 
tion, such as eccentricity, boundary curvature or the ex- 
istence of particular island chains, do not transparently 
determine emission characteristics in this regime. 




FIG. 1: Emission patterns are shown for two perturbations of 



a resonance of Vo(r ) = 



1) , with the generic perturba- 



tion ex in (a) and the integrable perturbation ex in (b). In 
each case the thick line shows the envelope predicted by first- 
order canonical perturbation theory and the thin line shows 
the angular dependence of |V>nm| 2 , normalised to have unit 
average, with n — 2 and m = 17 for e = —1/40 and h = 1/40. 
The respective Poincare plots, shown in (c) and (d), are de- 
fined by setting p r = 0. Even though the perturbations are of 
similar strength in each case, the primary islands near L — 
appear only at second order in the nonintegrable case and are 
much smaller, growing as e rather than \fe. 

We note that a "surprising" directionality in the emis- 
sion of weakly deformed resonators has previously been 
pointed out in Ref. [10] (see also [15, 16]), providing mo- 
tivation for the current study. However, existing analysis 
of such directionality has been based on the structure of 
real phase space and the role played, for example, by spe- 
cific island chains. This seems distinct from the mech- 
anism proposed here, which applies at deformations so 
small that real ray dynamics might hardly be affected. 
We will see that emission patterns can nevertheless be 
strongly modulated because there is an essential topo- 



2 



logical difference between the complexified dynamics of 
nonintegrable and integrable perturbations, which essen- 
tially amounts to a symmetry-breaking of the rays es- 
caping to asymptopia. These escaping rays are real in 
the integrable case but are slightly complex for noninte- 
grable perturbations and the resulting complex eikonal 
phase strongly modulates the outgoing wave. 

As a simple model of a deformed resonator we con- 
sider the Schrodinger equation for a perturbation V^(x) = 
V (r) + eV\ (x) of a central potential well V (r) . This 
is analogous to a scalar optical problem \7 2 ip(x.) + 
k 2 n 2 (x)tp(x) = with refractive index and wave num- 
ber defined by k 2 n 2 (x) = 2(E — V)/h 2 (assuming unit 
mass) . Of course, optical resonators have sharp material 
interfaces which should be modelled by step potentials in 
the Schrodinger equation. We present numerical evidence 
that a similar mechanism is at play for such systems in 
Fig. 2, where strong directionality is found in emission 
from perturbations of a circular resonator (Fig 2(a) and 
2(b)) which are so weak that there is little sign of nonin- 
tegrability in Poincare plots (Fig. (2(d)). We confine our 
detailed analytical work in this Letter however to smooth 
(or more accurately, analytic) potentials. The underlying 
idea can be set out more transparently this way and we 
will treat step potentials more properly in a future pub- 
lication. For similar reasons, we will present details for 
the two-dimensional case only and assert that it will be 
obvious how to adapt the discussion to three dimensions. 
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FIG. 2: Emission patterns are shown for two perturbations of 
a circular resonator with refractive index n = 2. An elliptical 
deformation is shown in (a) and a generic perturbation r = 
1+e cos 3x is shown in (b). In each case, e = 1/400, where the 
elliptical perturbation is r ~ 1 +e cos 2\ at first order and the 
resonance has m = 50 and Re(fc) as 33.19. A two-dimensional 
illustration and a Poincare-Birkhoff plot are shown for the 
generic perturbation in (c) and (d) respectively. In (c), the 
intensity of the emitted wave has been exaggerated to make 
it visible. 

We consider a Gamow-Siegert state tp nm (x), defined 
as a solution of the Schrodinger equation for a complex 
energy E = Eq — iF/2 which satisfies radiating bound- 
ary conditions at infinity. Here n and m denote radial 
and azimuthal quantum numbers of the resonance in the 
limit of zero perturbation. In the short-wavelength limit 



we can associate this state with an invariant torus of clas- 
sical rays moving in the plane configuration space in an 
annulus bounded by inner and outer caustics F% and T2, 
as illustrated in Fig. 3. 




FIG. 3: Real and complex rays are described in the complex 
angle plane as shown. Pi, T2 and T3 are caustics and Ci, 
C2 and C3 are contours in the complex 9 plane which define 
loops on these tori as shown at left. 

Immediately outside T2 is a band of complex rays 
where ip nm (x) decays exponentially in the radial direc- 
tion. We argue below that this band has an outer caustic 
T3 where it joins with a family of rays escaping to infin- 
ity (see Fig. 3) . It is these escaping rays which determine 
the emission pattern. The wavefunction outside T3 has 
an eikonal form 

^„ m (x) « A(*y s ^/ h (i) 

where S*(x) is the action function of the escaping rays. 
We expand the action in the form 

5(x) = S , (x)+e5i(x) + -.. , 

where Si (x) is determined explicitly below using canon- 
ical perturbation theory. For integrable problems the es- 
caping rays are necessarily real [17] and so therefore is 
VS*. This implies that any angular variation in the mag- 
nitude of the emitted wave derives from the amplitude 
A(x) and is slight if e is small. If the escaping rays are 
complex on the other hand, as we will argue they are in 
the generic nonintegrable case, then the emission pattern 
is dominated by the imaginary part of the action 

|i m (x)| 2 «e- 2dmS 'W/ fi 

at leading order. In this case a perturbation of the order 
of the small parameter h suffices to dramatically alter 
the envelope of the outgoing wave. We now set out the 
details of this calculation. 

We let (J, 9) and (L, cj>) respectively denote the action- 
angle pairs of variables for the radial and azimuthal de- 
grees of freedom in the unperturbed limit and let uj 
and f2 denote the corresponding frequencies. Because 
the actions are fixed by the torus quantisation rules 
I = (n + h)H and L = mh at leading order, we will sup- 
press the actions notationally in much of what follows, 
writing r{9) instead of r(9,I) for example. We denote 
by (r, x) the standard polar coordinates in the plane and 
point out that \ and </> are distinct but related by a shift 

X = + l(&), 
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where 

We adapt the convention that 9 = on the outer caustic 
T 2 of the real torus (Fig. 3). To characterise emission 
we must construct the complex rays in the classically 
forbidden region outside r 2 . We first describe how this 
is achieved in the unperturbed limit. 

Starting on the caustic r 2 where = 0, and letting 
move down the imaginary axis (the rays which define 
decaying WKB solutions are in the lower-half plane) , we 
find that the radial coordinate r{&) evolves periodically. 
Denote the imaginary period by 2iQ, so that 



where 



r(9 + 2i&) = r(0). 



(2) 



As angles (0, x) respectively range over the imaginary 
and real axes, a two-dimensional complex extension of the 
real torus is swept out which itself has the topology of a 
torus and on which radial momentum p r is imaginary and 
r is real. This second torus describes a band of evanescent 
decay immediately outside the resonator. It has an inner 
caustic r 2 , where it touches the real torus, and an outer 
caustic T^. The caustic corresponds to = — i<d. If 9 
evolves horizontally from — i<3 (contour C3 in Fig. 3) then 
we describe a family of real orbits escaping to infinity. 

The crucial feature in this description (of the unper- 
turbed case) is that the orbits escaping to infinity are 
real. This is related in the following way to the fact that 
the radial coordinate and other phase space functions are 
biperiodic functions of 6 — in addition to the usual real 
period r(9 + 2tt) = r(0) there is an imaginary period ex- 
pressed in (2). We have stated that the escaping rays are 
launched from the contour C 3 in the complex plane (see 
Fig. 3). They are real if they coincide with their complex 
conjugates, which are found along the conjugate contour 
C3. Since C£ = C3 + 2z6, the rays are therefore real if 
there is an imaginary period, as claimed. We will now 
show that if the rays are approximated using canonical 
perturbation theory, this second period is destroyed at 
first order and the rays therefore become complex. 

We use a type-two generating function -F 2 (0, <fi, I, L) = 
9I + 4>L + eG(9, <j>, J, L) to generate the transformation to 
action-angle variables (9, (j>, I, L) for the perturbed sys- 
tem. At first order, G(9, (j), T, L) provides the leading 
correction to the phase of Vw(x) and is a solution of 



dG BG 



(3) 



where U(9, <fr) is the oscillating part of V\ (x) , expressed in 
action-angle variables [14]. As in the unperturbed prob- 
lem, (I, L) are fixed by torus quantisation conditions and 
are suppressed notationally from now on. 

The perturbing potential U(0, (j)) inherits a biperiodic 
structure from the radial degree of freedom as follows. 
We first note that 

7(0 + =7(0) 



L 



r 2 (iu!T) 



dr 



is real. The potential is periodic with respect to (0, x) ~* 
(9 + 2iQ,x)i which, in terms of dynamical angles, gives 

U(9 + 2iQ, (j) + 2i$) = U{9, 0). 

Let the perturbing potential have the form 

U{r, x)=J2 U k(r)c tkx = J2 «k(*)e ifc * 

k k 

where u k {9) = U k (r(9))e lk "> (e l We substitute 

G(M)=£ffk(0)e <fc * 
k 

in (3) and solve the resulting equation for <?fc(0) to get 

-ikCl9 / uj 



9k(9) = c" 



ikne/o 



9k(0) 



to 



tkne ^u k (9')d9' 



The integration constant t/fc(0) is fixed by imposing the 
period g k {9 + 2n) = g k (9) on the real axis. The result- 
ing generating function is single-valued on the real axis, 
where it gives the real torus between caustics Ti and T 2 . 

Emission patterns are dominated by the imaginary 
part of G{0,<t>) along C 3 . In light of C* = C 3 + 2i6, 
this can be written ImG(0, <f>) = AG(6,(p)/(2i), where 

AG(0, <jy) = G(0, <jy) - G(9 + 2iO, + 2i$). 

It can be shown that 

AG(9,0) =i KQ + iK(9,<j)), 

where 

r 2&/w 



K 



— / M O (0)d0 = / Uq{iljJt)At 
Jo Jo 



is a real constant and 

kjtO C 

The contour C is a clockwise circuit of the rectangle at 
right in Fig. 3 (a derivation of an analogous result in the 
context of double- well splittings can be found in [14]). 
It is convenient to express the result using the angles 

a(6,cp) = (j) + i^--(9 + ie), 
to 

m = ^(0 + ,e ) + 7(0 + ie) = i£ Q ^d0'. 

Note that a is real and constant along the escaping rays 
of the unperturbed problem and is the value of the polar 
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angle % at which a given ray touches the caustic T3. We 
also have \ = a + Pi so /3 measures the angular deflection 
of the escaping ray as it moves away from T 3 . Then 

fe#0 

where 

The main conclusion of this Letter is that the asymptotic 
intensity of a WKB mode (1) is modulated by a function 

M (a) = e- rf (")/ R (4) 

which can vary strongly with \ even f° r e — 0(h). 

We illustrate this result numerically for the potential 
Vo(r) = — (r 2 — l) 2 . Although this potential does not van- 
ish at infinity, it defines Gamow-Siegert states entirely 
analogous to those of asymptotically free problems and 
is used because solutions are easily obtained for it using 
complex rotation in a harmonic-oscillator basis. None of 
the calculations above are substantially affected by this 
choice. We consider two perturbations, Vi(x) = x, which 
behaves generically, and Vi(x) = x 2 , for which the total 
potential can be shown to be separable in elliptic coor- 
dinates and which therefore provides us with an exam- 
ple of an integrable perturbation. The respective angu- 
lar distributions of the outgoing intensity, normalised to 
have unit average, are shown in Figs. 1(a) and 1(b). The 
heavy curves show the modulation function (4). Note 
that for the integrable perturbation we necessarily find 
that K(a) = and the small modulation in the numerical 



solution is a result of corrections to the amplitude A(x), 
which have not been taken into account. In the nonin- 
tcgrable case, the modulation argument can be written 
K(a) = 2ki cos a and provides a good description of the 
numerical emission pattern (Fig. 1(a)). 

A similar modulation is seen in the numerical re- 
sults for a deformed cavity shown in Fig. 2. The 
dielectric ellipse is a special case — it is not sep- 
arable but the real and complex rays can nonethe- 
less be calculated analytically. The envelope shown is 
M(a) — exp[2efc(n 2 — 1) n 2 p 2 — 1 cos 2a] , where p = 
m/k, whose derivation will be described in a future pub- 
lication. The envelope for the generically nonintegrable 
perturbation in (b) has not been calculated analytically 
but has been found numerically to be well described by a 
function of the form A4(a) = cxp[eka(p) cos3a], with a 
typical case shown in the figure. This is consistent with 
the mechanism set out in this Letter. 

In conclusion, canonical perturbation theory married 
with complex WKB approximation can successfully de- 
scribe strongly directional emission from weakly de- 
formed resonators. The underlying mechanism is not 
transparently related to the structure of real phase space 
and there is a stark qualitative difference between inte- 
grable and nonintegrable systems, even when the signa- 
tures of nonintegrablilty are slight in real phase space 
[18]. Despite this sensitivity to dynamical detail, simple 
analytical formulas can be given for the observed output. 
The theory has been presented for smooth potentials but 
a similar mechanism is expected to govern sharp cav- 
ity problems and should provide a useful analytical tool 
to understand the evanescent field outside optical res- 
onators, which is important to applications such as lasers 
and sensors [2]. 
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